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ABSTRACT 

In the redshift range z = 0—1, the gamma ray burst (GRB) redshift distribution 
should increase rapidly because of increasing differential volume sizes and strong evo- 
lution in the star formation rate. This feature is not observed in the Swift redshift 
distribution and to account for this discrepancy, a dominant bias, independent of the 
Swift sensitivity, is required. Furthermore, despite rapid localization, about 50% of 
Swift and pie-Swift GRBs do not have an observed optical afterglow and 60-70% of 
GRBs are lacking redshifts. We employ a heuristic technique to extract this redshift 
bias using 69 GRBs localized by Swift with redshifts determined from absorption or 
emission spectroscopy. For the Swift and HETE+BeppoSAX redshift distributions, 
the best model fit to the bias at z < 1 implies that if GRB rate evolution follows the 
SFR, the bias cancels this rate increase. We find that the same bias is affecting both 
Swift and HETE+BeppoSAX measurements similarly at z < 1. Using a bias model 
constrained at a 98% KS probability, we find that 72% of GRBs at z < 2 will not have 
measurable redshifts and about 55% at z > 2. To achieve this high KS probability 
requires increasing the GRB rate density at small z compared to the high-z rate. This 
provides further evidence for a low-luminosity population of GRBs that are observed 
in only a small volume because of their faintness. 
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1 INTRODUCTION 

Because GRBs are the brightest transient cosmological 
events measurable, they are a unique probe to the high- 
z Universe and provide an opportunity to study their lo- 
cal interstellar environment, host galaxies and the inter- 
galactic medium. Multi- wavelength observations of GRBs 
have confirmed that a significant fraction of long GRBs Q 
are associated with the collapse of sho rt-lived massive stars 
|Hiorth et al.ll2003l ; IStanek et alfeOOof ). This provides strong 
motivation to use the measured GRB redshift distribution as 
a tool to constrain the star formation rate (SFR) and GRB 
rate evolution. However, it is clear that the relation between 
the evolving SFR and the GRB redshift distribution is more 
complex than originally assumed. 

Studies of the host galaxies of supernovae reveal that 
they are more probable in brighter, massive spirals, while 
GRBs are more typ ically in smaller irregular host types 
ijFruchter et al-lfeoOrj ). A possible explanation for the dif- 
ference in host galaxy types is a preference for GRBs to 
form f rom relatively lower-m etallicity core-collapse super- 
novae (lYoon fc Langerll2005l ). This argument is compatible 
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with the 'collapsar' model l|Wooslevlll993l ), which postulates 
that massive GRB progenitors may be preferred in lower- 
metallicity environments. The inferred metal abundances in 
the GRB host galaxies are on average smaller than those 
found in ultra-luminous dusty galaxies (e.g. Le Floc'h et al. 
2006). Although the SFR is lower in the GRB host galaxies, 
it still may comprise about 2/3 of the global SFR. Hence, the 
observed GRB redshift distribution is possibly tracking the 
SFR in the faint end of the galaxy luminosity distribution. 

Up to 2007 November, more than 80 GRB redshifts have 
been measured from rapid localization by the Swift satellite 
and follow-up spectroscopy of the optical/NIR afterglow and 
host galaxies by large ground-based telescopes. Optical/NIR 
afterglows have been found for nearly 60% of Swift localized 
GRBs but only about 30% of Swift GRBs have measured 
redshifts, implying a very incomplete sample. The reasons 
why a significant fraction of localized GRBs have no de- 
tectable afterglow is an ongoing controversy. It is known 
that GRB afterglow radiation is sub ject to several a bsorp- 
tion processes in the optical band I Gou et ail [2004) . The 
presence of dust and gas, though, does not obscure gamma 
rays, irrespective of whether this materi al is within the host 
galax ies or in the intergalactic medium l|Blain fc NataraianI 
|2000|) . 
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2 OPTICALLY DARK GRBS 

GRB optical afterglows may be dark due to obscuring dust, 
very high redshifts or rapidly decaying transients, a phe- 
nomena which is commonly referred to as a dark burst-see 
iMirabel fc Rodriquesl (2003) and references therein. A sub- 
class of 'in trinsically dark' GRB afterglows has been argued 
against by M cszaros et alJ (|2005h . They show that selection 
effects are strong enough to explain the proportion of GRBs 
deemed to be optically dark. Such 'dark' bursts may be a re- 
sult of dust extinction in the local ISM of the host along the 
line of sight. Evidence for this is shown in a recent study of 
the host galaxy of GRB 030115 l|Levan et al.ll200rj ). Its op- 
tical afterglow was fainter than many upper limits for other 
bursts, suggesting that without early NIR observations it 
would have been classified as a dark burst. Both the colour 
and optical magnitude of the afterglow are probably affected 
by dust extinction and indicate that at least some optical af- 
terglows are very faint owing to dust along the line of sight. 
Furthermore, the GRB host galaxies are on average smaller 
and fainter than supernova host galaxies, so it is more dim- 
cult to identify many GRB hosts. 

In another study. iRoming et al.l (|2006l ). used very early 
observations of Swift GRBs to investigate the most proba- 
ble cause of optically dark GRBs. They find that ~ 25% of 
the bursts in their sample are extincted by Galactic dust, 
~ 25% are obscured by absorption in the immediate GRB 
environment and ~ 30% are most likely lost by Ly-a blan- 
keting and absorption at high redshift. So it is highly likely 
that the dark bursts result from a combination of extinc- 
tion factors relati ng to the GRB envi ronment, host galaxy 
type and redshift. ISchadv et al.l (|2007l ) estimate that GRBs 
without an afterglow detected in A > 5500 Angstrom have 
average optical extinctions eight times those observed in the 
optically bright population of GRBs. They point out that 
extinction could account for the non-detection of GRB af- 
terglows by the UVOT onboard Swift. 

Given the incompleteness of the Swift GRB redshift 
sample, one strategy for using the current redshift distribu- 
tion as a high-z probe is to remove bursts from the sam- 
ple where optical afte rglo w observation was sub-opt imal . 
iJakobsson et al.l (120061 ) and lTanvir fc Jakobssonl (|2007h used 
a set of criteria to select a sub-sample of GRBs based on 
Galactic foreground extinction and Sun and Moon positions 
relative to the burst. Despite this culling they find that the 
redshift completeness of Swift GRBs could be increased only 
to about 50-60%, implying that a combination of other se- 
lection effects are severely affecting the sample. 

In this study, we use a different approach, that focusses 
on investigating how the redshift sample is modified by a 
combination of unkown selection biases. We show how a 
global selection function can be extracted from GRB red- 
shift data selected with minimal assumptions. We do not 
attempt to model the selection function based on strong 
physical arguments, because of the large uncertainty in the 
multiple biases that comprise the selection function. Regard- 
less of this limitation, the form of the selection function ex- 
tracted from our analysis provides a means to estimate how 
the dominant biases may be modifying the redshift distribu- 
tion. To achieve this aim we employ a subset comprised of 
69 Swift redshifts obtained from well-localized GRBs with 
Tgo > 2s and accurate spectroscopic redshifts. We have not 



included photometric redshifts or upper and lower limits 
for redshifts. The sample includes the controversial GRBs 
060614 (z = 0.12) and 060505 (z = 0.089), which were long 
duration events without accompanying supernovae. 

The goal of this study is twofold: firstly to constrain a 
model for the combined biases using a heuristic approach 
with minimal assumptions. Nonetheless, we do include a 
very crude model for the so-called 'redshift desert' at z « 
1 — 2 - a region where it is difficult to measure redshifts be- 
cause of the lack of obse rvable strong emission or absorption 
lines (|Fiore et al.ll2007h . We point out that this bias is re- 
stricted to those GRBs with optical afterglows-i.e. 50-60% 
of the total Swift GRB sample, implying a combination of 
other astrophysical and observational selection effects mod- 
ify the sample. To account for the remaining missing fraction 
of redshifts, we employ a parametrized function that scales 
the fraction of GRB afterglow detections at z = 0.1 — 6. 

Secondly, we will examine if the fractional contributions 
of individual selection biases (e.g. Roming et al. 2006) are 
reconcilable with the global selection function extracted in 
this study. The elegance of this technique is that the selec- 
tion function tracks the evolution of the dominant biases in 
different redshift ranges. The selection function extracted in 
our analysis is constrained using a KS test applied to the 
chosen subset of GRB reshifts and a bias-corrected proba- 
bility distribution model. This is a first step toward using 
the distribution of 'missing' GRB redshifts as a potential 
new probe of the evolution of GRB environments. 



3 THE GRB REDSHIFT SELECTION 
FUNCTION 

An inspection of the Swift redshift distribution provides an 
insight into the form of such a function. Firstly, in the red- 
shift range 0—1, the distribution should increase rapidly be- 
cause of increasing differential volume sizes. Assuming that 
GRB rate evolution follow the global SFR in a 1:1 manner, 
recent observation-based models of the SFR at z < 1 imply 
a rapid increase in the GRB redshift distribution in the same 
range. This feature is absent from the present distribution. 
Secondly, the distribution shows a distinctive drop in num- 
bers at z ~ 1 — 2, and an increase at z > 2. Figure 1 plots 
the Swift normalized redshift distribution (shaded region), 
along with a model distribution. 

To model a complete GRB redshift selection function 
requires knowledge of the component biases. One is the sen- 
sitivity of the satellite detector assuming a distribution of 
source luminosities. For Swift, there are also complex de- 
tector triggering algorithms that may introduce a bias, but 
these effects have not been modelled completely and we do 
not attempt to account for them here. Since redshift mea- 
surement depends on obtaining optical emission or absorp- 
tion spectra, it is clear that any selection of optical after- 
glows will affect the redshift distribution. Because about 
50% of GRBs are optically dark, there is a dominant selec- 
tion effect operating on the optical afterglow distribution, 
as discussed above. 

Of the GRBs that do have observed optical afterglows, 
only about 60% have measured spectroscopic redshifts. It 
is apparent that other biases related to to the acquisition 
of spectroscopic redshifts may play a part. The 'redshift 
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desert', z ~ 1 — 2, is a region where it is difficult to mea- 
sure redshifts because of the lack of strong emi ssion or ab- 
sorpt ion lines covered by optical spectrometers l)Fiore et al.l 
120071 ). Strong emission lines such as Ha, H/3, OIII and Oil go 
out of the observed range at z ~ 1.1 while Lyman-Q enters 
the range at z ~ 2.1. The strong absorption feature of the 
Mgll (A = 2796, 2803) doublet is contaminated by telluric 
features at z > 1.5. 

The magnitudes of the individual biases from these ef- 
fects are expected to be very different. For instance, dust 
extinction may well be a dominant effect given that a signif- 
icant fraction of GRBs are optically dark. Our approach is to 
estimate a 'global' selection function by weighting the red- 
shift sample by the flux-limited sensitivity of the detector 
and volume-number dependence from the redshift sample. 
The geometric component is well known but the flux-limited 
sensitivity is less so because of the uncertainty in the GRB 
luminosity function. 

For self consistency we assume a single luminosity 
function (LF) that encompasses the high luminosity (HL) 
GRBs and the so-called anomalous low-luminosity (LL) 
GRBs, those with luminosities < 10 49 ergs s _1 , detected by 
HE TE+ B eppoSA X and Swift. We adopt the d o uble power- 
law LF, 4>(L), from lGuetta, Piran fc Waxmanl ([20051) . with 
slopes —0.1 and —2, with luminosity limits that encompass 
the extremely LL GRB 980405: (L max ,L min ) = (3xl0 53 , 3x 
10 46 ) ergs s _1 . The fraction of detectable GRBs, or flux- 
limited selection function (those observed with peak flux 
> Fi im ), is 

i/>fiux(z) = / 0(L iso )dL iso , (1) 

J £lim(flim,z) 

where Li so is an isotropic equivalent luminosity and 
£iim(-Fiim, z) is obtained by solving Fiim(Ln m ,z) = 
L\i m / (47tdi,(z)), with di J (z) denoting the luminosity distance 
of the burst. The flux sensitivity limit, -Fn m , of the detector 
BAT onboard Swift is about 1 x 10 -8 ergs cm -2 s _1 and for 
a combined HETE+BeppoSAX sensitivity we take a mea n 
of 4 x 10~ 8 ergs cm -2 s -1 , approximated from lBandl (2006). 



3.1 The GRB redshift probability distribution 
function 

The main redshift distribution biases are represented as 
the product i/'fl ux (z)'i/'redahift(«) of two dimensionless selec- 
tion functions, where ?/> re dshift (z) describes an evolving bias 
comprised of both spectroscopic and astrophysical compo- 
nents; ^redshift (z) is uncertain, but plays a fundamental role 
in shaping the red shift probability distribution function — 
sec Oowar d (2007) for a more detailed derivation. The GRB 
probability distribution function is expressed as: 



P(z) =N 1 ipa ux (z)i>rcdshiit(z)e(z) dV J'^ ^ ^ 
with normalizing factor 



N = 



V'flux(«)V'red s hift(z)e(2:) 



dV(z) dz 
dz (1 + z) 



(2) 



(3) 



The (l + z) -1 accounts for cosmological time dilation, dV(z) 
is the volume element in a flat-A cosmology (Qm = 0.3, 
SIa = 0.7) with Ho = 70 km s _1 Mpc -1 and z max is the 
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Figure 1. Plot of the normalized distribution of 69 spectroscop- 
ically measured redshifts obtained following Swift localizations, 
shaded, and probability distribution model for the same distri- 
bution, with no correction for biases. The m odel assumes a GRB 
rate e volution that follows the SFR model of lHopkins fc Beacoml 
(2006). It peaks at z = 1 — 2, where the distribution tracks the ge- 
ometric increase in numbers. At z > 2, the rapid fall-off is caused 
by the reduction in flux sensitivity of Swift and decreasing cos- 
mological volume shells. 



maximum observable redshift. We assume that the dimen- 
sionless GRB rate evolution factor, e(z), follows the evolving 
SF R density, p,(z). 

iHopk ms fc Beacoml (|2006T ) have constrained the SFR 
by fitting to ultraviolet and far-infrared data obtained up 
to 2006. When these data are taken together, they show 
that it is possible to constrain the SFR out to z ~ 6, with 
especially tight constraints for z < 1. Equation ((4]) shows 
the functional form they employed in their fit to the SFR: 



e(z) 



P.(*) 

MO) 



1 + a\z 

+ (z/a 2 ) a 3 



(4) 



where the evolution function e(z) is defined so that e(0) = 1. 
We use the best fit parameters for e(z) based on this SFR 
model: (or, a 2 ,a 3 ) = (7.6,0.0059,4.3). 

Combining (1 + z)~ 1 dV/dz with the flux- limited satel- 
lite selection bias, ^/>flux(z), yields a resulting distribution 
that is strongly peaked at z = 1 — 2. Figure 1 plots the 
distribution described by equation (J5J) assuming t/Wshift (z) 
is constant (i.e. no bias correction). The Swift GRB red- 
shift distribution, comprised of 69 spectroscopic redshifts, is 
clearly different in shape, highlighting a redshift dependent 
selection bias other than the flux sensitivity of the satel- 
lite. The most obvious features are the lack of numbers at 
2 ~ 1 — 2 and the relatively high numbers at z > 3 compared 
to * = - 1. 



4 ANALYSIS TECHNIQUE 

Because the probability distribution function @ is sepa- 
rated, the combination of flux-limited selection bias and 
geometric components, ipa ux dV/dz, can be used to weight 
the observed redshift distribution. Furthermore, because the 
SFR is reasonably constrained at z = — 1, and assuming 
that e(z) follows the SFR, this can also be used to weight 
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the observed redshift distribution. If V>redshift(z) evolves with 
redshift, the residuals that remain after weighting the ob- 
served distribution will not be uniformly distributed. The 
aim of the analysis is to construct models for Vredshift (2) 
that minimize the variance of the residuals after weighting 
the observed redshift distribution. As a further constraint, 
we calculate cumulative distributions for the observed and 
bias-corrected model probability distributions and apply a 
KS test to quantify their compatibility. 

A further observational constraint on Redshift (z) is the 
fraction of measured redshifts out of the total GRBs local- 
ized by Swift, namely about 30%. We approximate this func- 
tion for the redshift detection efficiency so that it scales ac- 
cording to the fraction of GRBs with measured redshifts. 
The redshifts are split into 3 ranges for the Swift and 
HETE+BeppoSAX samples to help identify the effect of dif- 
ferent redshift biases. The first, z — — 0.1, is dominated 
by LL GRBs that occur at a possibly higher rate than the 
HL GRBs. In the following analysis we constrain the relative 
rates of the LL bursts to the HL high-z population in this 
small redshift range. 

The rapid decrease in redshift numbers at z < 1 from 
the observed distribution is similar to the inverse of the slope 
of e(z) in the same range. It is difficult to determine if there 
is a direct relationship between this bias and the SFR. Dust 
obscuration may play an important role in modifying the 
afterglow detection rate, but given the substantial uncer- 
tainty in the links between dust obscuration, SFR evolution 
and GRB redshift detection, we take a heuristic approach. 

We include a crude model for the redshift desert. From 
figure 3 of Bloom (2003), the probability of obtaining promi- 
nent emission and absorption lines at 1.5 < z < 2 is reduced 
by at least 50%. The same figure shows that at 2 < z < 3.5 
the probability of detection increases significantly. We em- 
ploy a parametrized toy model for the redshift desert that 
scales the probability of detection at 1 < z < 2. It is essen- 
tially unity outside of this redshift range and decreases to a 
minimum at z = 1.6: 



ip d (z; A, B) = 1 - (1 - A)sech 2 [B(z - 1.6] , 



(5) 



where A is the minimum value and B defines the width of 
the desert. Figure [2] plots ip d (z) with (A, B) = (0.15,3.00). 
In this model, redshift determination is affected only in 1 < 
z < 2, where the desert is prominent. 

We separate the contributions of ^redshift (z) into two 
parts. The first, comprised of LL GRBs, allows for a rela- 
tively higher rate density than the HL GRBs. At z < 0.1, the 
Rll parameter scales the GRB rate density of LL GRBs, 
all observed at small redshift, to the HL high-z bursts: 



</>ll(z) = Rll for z < 0.1 



(6) 



For z > 0.1, Vwshift (z) is constrained by the Swift HL GRB 
redshift distribution, and assumes the form ip d (z)/e(z). Two 
broad redshift ranges are defined, the first one is dominated 
by tpA^z) and l/e(z), the second by l/e(z) and a redshift 
dependent power law of the form z 1 : 



1pHh(z) 



aipd(z) I 'e(z) 
Pz"<Mz)/<z) 



for < z < 1.85 
z > 1.85 



(7) 



where z — 1.85 is chosen to ensure that the model joins 
smoothly at the boundary between the two redshift ranges. 
The parameters a, (3 and 7 are constrained by finding the 




Figure 2. Plot of the the toy model for the 'redshift desert', equa- 
tion J5}, assuming parameters (A, B) = (0.15,3.00) and centred 
on z = 1.6. The effect of this suppression on redshift determina- 
tion is negligible outside of z = 1 — 2. 



highest KS probability using the weighting technique de- 
scribed above. The total bias function is defined as: 



^redshift (2) = ^LL(z) + ^HL (z) 



(8) 



We apply the same technique and parametrized 
form to the sample of 33 GRB redshifts localized by 
HETE+BeppoSA X, assuming the m ean flux threshold of 
these detectors. iFiore et al.l (|2007l ) point out that the 
median delay time of optical follow-up for Swift and 
HETE+BeppoSAX GRBs is 15 minutes and 3.5 - 14 hours 
respectively. This results in a smaller proportion of faint af- 
terglows triggered from HETE+BeppoSAX, so that a bias is 
introduced against measuring redshifts from high-z GRBs - 
i.e. those with faint optical afterglows. To account for this 
bias, we employ a power law of the form Pt = z~ s at z > 1, 
and constrain 8 using the KS test. 

Given that i/> ro dshift can be constrained, the fraction of 
HL GRB redshifts that are not observed as a function of 
redshift, Fdark(z), can be calculated from 



Fdark(z) = 1 — Z~ 



ipHh(z')dz' for z ^ 5 



(9) 



We assume that Fd ar k(z) represents the contributions 
from optically dark GRBs and those with observed after- 
glows, but with no measured redshift. From the Swift sam- 
ple we have selected, the total fraction comprised of both 
contributions equates to 60 — 70%. It is interesting to con- 
sider the affect on Fdark{z) of using a more selective sample. 
For example, a sample excluding any bursts located near the 
Sun or Moon increases the completeness of the GRB redshift 
sample. The effect of fractional completeness on Fdark(z) is 
to scale linearly the function but not change its shape unless 
redshifts are omitted from the sample. Because the redshift 
sample we employ is statistically small, we have chosen not 
to reduce the sample size further. In fact, doing so would be 
contrary to our main goal, which is to study how the biases 
modify the redshift distribution. 
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5 RESULTS 

We apply the weighting technique to the Swift redshift dis- 
tribution using the parametrized model above to constrain 
the fitting parameters in the redshift ranges 2 = 0.1 — 1.85 
and z > 1.85. At the 90% KS probability level, the free 
parameters are constrained to (Rll, oi, /3, 7) = (20-90, 1.3- 
1.8, 0.35-0.45, 1.45-1.55). Cumulative distributions of the 
Swift redshift distribution and a bias correction model are 
plotted in figure [3] and [4] We find the highest KS probabil- 
ity of 98% using (Rll, a, /3, 7) = (40, 1.7, 0.38, 1.5). For 
the HETE+BeppoSAX sample we find similar values for the 
above parameters, except for Rll ~ 30, slightly smaller 
than that obtained using the Swift distribution-see figure 4. 
To be compatible with the HETE+BeppoSAX distribution, 
the bias model has to include the additional suppression of 
detections at z > 1 with S — 0.4. 

In the volume bounded by z ~ 0.1 for the Swift and 
HETE+BeppoSAX distributions, the optimal KS probabili- 
ties imply a rate increase of 20 — 90 times that of the higher- 
z HL GRBs. This could be evidence for a separate GRB 
population with a luminosity distribution biased towards 
mu ch lower l umin osities, compared to the HL population- 
see [Oowar 1 (120051) for an early a ttempt at modelling this 
population and Liang et alj {2007) for more detailed analy- 
sis. Alternatively, the LL GRBs may represent a component 
of a single population characterized by nearly isotropic emis- 
sions. 

Here, we assume a single LF encompasses all GRB lu- 
minosities. The relative increase of the GRB rate in small 
z may be attributed to two factors: a larger mean GRB jet 
opening angle, 9, and the small volume that the GRB de- 
tectors are sensitive to. The beaming factor, HL (0) = 
1— cos 9, defines the number of GRBs not observed for ev- 
ery one that is observed. We can estimate 9 for the nearby 
GRBs via f~l{9) = & l (6)/Rll, where f£(0) and f HL {9) 
are the beaming factors for the LL and HL GRBs respec- 
tively. Using Rll ~ 40 and an estimate of 10° for th e HL 
GRB beaming angle l|Guetta. Piran fe Waxmanl I2005T ) . we 
obtain 9ll ~ 70°. This is compatible with observations of 
the late-time radio afterglow of GRB/XRF 060218, with a 
j et opening angle that appeared to be much less collimated 
{Soderberg et al.ll2006T ). 

Using the bias correction model ^redshift (2) with (a, 0, 
7) = (1.7, 0.38, 1.5) and equation ©, we calculate the frac- 
tion of unobserved redshifts for GRBs localized by Swift. We 
find that 72% of GRBs redshifts are missing at z < 2, 55% 
at z > 2 and a total fraction « 60% over all redshift. This 
result is independent of the physical origin of the biases. 

Figure [5] reveals some interesting trends in optical af- 
terglow and redshift non-detection, as a function of redshift. 
It is clear that the strongest evolution of optical afterglow 
extinction is occurring at z — — 1, outside of the redshift 
desert at 2 ~ 1.5 — 2. For 2 > 2, the effect of the bias 
is reducing with redshift, implying that the probability of 
measuring a redshift is increasing with redshift. 



6 DISCUSSION 

We show that as the total rate of GRBs increases with the 
SFR at 2 = — 1, the observed fraction of GRB redshifts 
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Figure 3. Top panel — plot of the cumulative distribution of Swift 
GRB redshifts and a model distribution not corrected for biases. 
Lower panel — the best fit model that includes bias correction 
using parameter values for Vredshift(z) °f {Rll, &, f3, 7) = (40, 
1.7, 0.38, 1.5). The corresponding KS probabilities are stated. 



decreases at a similar rate. This is a surprising result, and in 
effect implies that GRB rate evolution, assuming it follows 
the SFR, cancels from the distribution model. The fitting 
results for 2 < 1 imply that if GRB rate evolution is linked 
to the SFR, there must be a bias that compensates for this 
rate increase. One possibility is that GRB redshift determi- 
nation becomes more difficult in redshift ranges where star 
formation is prolific. It is difficult at present to untangle the 
relationship between the SFR and dust extinction of GRB 
afterglows. But we do know that the SFR evolves strongly 
at 2 = — 1, the range where most of the GRB redshifts 
are missing. Using the most probable model for VWshift(.z), 
Fdark(z) yields a total fraction of missing GRB redshifts 
~ 63%, a result supported independently by the observed ra- 
tio of GRB redshifts to GRBs of about 70%. Given that 55% 
of GRB optical afterglows are exting uished by host galax y 
dust and absorption at high redshift l|Roming et alj 2006), 
our results are suggestive of an extinction function that may 
compensate for GRB rate evolution tracking the SFR. 

An inspection of figure [5] at high 2 presents a puzzle. 
If the dominant biases that modify the redshift distribution 
are distance dependant, such as dust extinction, then why 
is the probability of observing an afterglow and measuring 
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Figure 4. Same as figure 3, but using the HETE+BeppoSAX 
redshift distribution with the same optimized parameter values 
shown above. To account for the relatively longer average time 
delay in optical afterglow follow-up compared to Swift, we include 
a suppression factor at z > 1 of z~ s . The optimal KS probability 
constrains 5 to 0.4. 

a redshift increasing with redshift. Several speculative ideas 
that are compatible with this result include: 

(i) GRBs are increasing in luminosity at high redshifts, 
resulting in a bias for selecting optical afterglows at z > 2. 

(ii) the GRB rate at z > 2 is evolving faster than the 
assumed global SFR in the same redshift range. 

(iii) an evolution in the type of dust-i.e. evolving opacity 
at high redshift. 

The first suggestion has been proposed before to explain the 
relatively large number of highly luminous GRBs observed 
at high z. We note that if we assume that the GRB rate does 
not follow the global SFR, but increases more gradually, the 
strong bias at z = 0.1 — 1 is reduced. This is obvious from our 
analysis, since GRB rate evolution that follows the global 
SFR is essentially cancelled out by the selection function. 
But we cannot rule out that the possibility that the global 
SFR model used in this study underestimates the 'true' rate 
at high z. 

Finally, our analysis supports other work finding in- 
creased GRB numbers at small z, corresponding to LL 
GRBs. For both the Swift and HETE+BeppoSAX distribu- 
tions we find a similar increased rate of GRBs is required at 
small redshift compared with high- 2 bursts. If the enhanced 
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Figure 5. Plots of the most probable redshift detection bias func- 
tions for GRBs localized by Swift and HETE+BeppoSAX using 
^redshift with (R LL ,a,/3,-y) = (40,1.7,0.38,1.5). The normal- 
ized bias function is the same for Swift and HETE+BeppoSAX 
apart from the power law with 8 = 0.4 at z > 1 for the lat- 
ter. Data, squares for Swift and triangles for HETE+BeppoSAX, 
are obtained from weighting the distributions by ipfl ux (z)e(z)(l + 
z)~ 1 dV/dz. The models include an enhanced rate density of 40 
at z < 0.1 to achieve a maximum KS probability > 98%. 

rate results from a nearly isotropic em ission component, the 
rates of HL and LL GRBs are similar l|Guetta fc Delia Valle 
l2007| y 

The challenge of understanding the dominant bias from 
the GRB redshift sample can be tackled in two ways. Firstly, 
small number statistics result in significant uncertainty in 
the parameters that comprise Vredshift(z). As the number 
of Swift redshifts approaches one hundred, the shape of 
V'redshift (z) at 2 = — 1 should be further constrained to 
the point where a formal fit to the data becomes more use- 
ful. Secondly, to achieve this goal requires constructing mod- 
els for 7/i r edshift (z) that relate to the physical processes that 
cause the bias. A significant improvement to the heuristic 
approach used here would be to implement a universal ex- 
tinction function for GRB afterglows. This would provide a 
test of whether such a model is reconcilable with the evolu- 
tion of biases reported here. 

It is especially critical to untangle the relationship be- 
tween dust extinction of the optical afterglow, GRB redshift 
determination, GRB rate evolution and the SFR. Is the re- 
duction in redshift numbers observed at z — — 1 somehow 
related to the SFR, or does this suppression originate from 
an entirely different class of biases? An ironic possibility that 
this work highlights is that the population of GRBs with- 
out measured redshifts may tell us just as much about GRB 
environments as the population with redshifts. 
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